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It is shown that the continuum level density (CLD) at unbound energies can be calculated 
with the complex scaling method (CSM), in which the energy spectra of bound states, 
resonances and continuum states are obtained in terms of L 2 basis functions. In this method, 
00 ■ the extended completeness relation is applied to the calculation of the Green functions, and 

the continuum-state part is approximately expressed in terms of discretized complex scaled 
continuum solutions. The obtained result is compared with the CLD calculated exactly 
£f} . from the scattering phase shift. The discretization in the CSM is shown to give a very good 

^ ' description of continuum states. We discuss how the scattering phase shifts can inversely be 

' calculated from the discretized CLD using a basis function technique in the CSM. 
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Recently, there has been much interest in nuclear structures of unstable nuclei, 
in which exotic nuclear structures have been revealed through the development of 
radioactive nuclear beam experiments. 1 - 1 It has been shown that for such nuclei, for 
example the so-called neutron halo nuclei, there are extremely weak binding ground 
states, and most of the excited states are in the continuum energy region. Therefore, 
to understand the exotic structures and excitations of these nuclei, it is necessary to 
study continuum and resonant states in unbound energy regions. 

The continuum level density (CLD) is expected to play an important role in 
relating experimental data and theoretical models for unbound states. Recently, 
Kruppa and Arai 2 ^ 4 ^ proposed an interesting method to calculate the CLD and 
argued that resonance parameters can be determined from the CLD. They start 
their investigation from the definition of the CLD, 5 ) 



A(E) = --Im[Tv [G(E)-G (E)}}, (M) 

7T 



where the full and free Green functions are given by G{E) = (E — H)^ 1 and 
Gq(E) = (E — Ho) -1 , respectively. Because the Hamiltonian H includes finite range 
interactions in addition to the asymptotic Hamiltonian Hq, the CLD expresses the 
effect from the interactions. When the eigenvalues (ej and e J , respectively) of H and 
Hq are obtained approximately within a framework including a finite number (TV) 
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of basis functions, the following discrete level density is denned: 

N N 

A N (E) = - e ^ " E ^ ~ e o)- (!- 2 ) 

i j 

Kruppa 2 ) employed a smoothing technique defined by the Strutinsky procedure 6 ) to 
calculate the continuous CLD, A(E), from its discrete form, A^(E). However, as 
discussed in Ref. 4) their results for the CLD exhibit a strong dependence on the 
smoothing parameters. We desire a more effective method to smooth the discrete 
quantities, or to discretize the continuum states. 

In this paper, we study a more direct method to calculate the CLD with no 
smoothing technique in the framework of complex scaling, 7 ) in which a basis function 
method is used to obtain not only bound states but also resonance and continuum 
states. The idea for the present method is taken from the extended completeness 
relation, 8 ) originally proposed by Berggren, 9 ) for bound, resonance and continuum 
states in the complex scaling method (CSM). Exact proofs of this extended com- 
pleteness relation for the CSM were recently given for a coupled channel system 10 ) 
and a single channel system. 11 - 1 Green functions can be expressed by using the ex- 
tended completeness relation in terms of discrete eigenvalues of the CSM with a finite 
number of basis functions. Because the complex scaled Hamiltonians H e and H® 
have complex eigenvalues, singularities, like the ^-function contained in Eq. (|l-2[) . 
are avoided and replaced by Lorentzian functions. Therefore, no smoothing pro- 
cess is needed. Furthermore, it is shown that A(E) for the CLD can be calculated 
independently of the scaling parameters in the CSM. 

Kruppa and Arai 2 )~ 4 ) applied the CLD to search for resonance parameters. How- 
ever, although parameter values for narrow resonances can be obtained using any 
method of CSM and CLD, it is not easy to extract parameter values for broad reso- 
nances with the CLD. Rather than obtaining such resonance parameters, it is more 
important to calculate the phase shift and/or S-matrix for the scattering states from 
the CLD. The CLD A(E) is related to the scattering phase shift 5(E) as 5 )' 12 ) 

, . 1 dS(E) , x 

4 < E > = ?-ET (1 ' 3> 

Therefore, once we confirm that the CLD calculated in the present method is con- 
sistent with A(E) obtained from the phase shift of the scattering solution, we can 
inversely calculate the phase shift by integrating the CLD obtained as a function of 
the energy from the eigenvalues of the complex scaled Hamiltonians. This implies 
that the phase shift can be calculated from discrete eigenvalues using a basis function 
method. 

We here demonstrate the reliability of this method by applying it to several 
two-body systems, including 4 He+n and a+a, which were studied by Arai and 
Kruppa. 4 ) Comparing the calculated continuum level density and the phase shifts 
in the CSM with the results obtained from the exactly calculated phase shifts and 
their derivatives (the CLD), we show that the extended completeness relation in the 
CSM is effective, and also that the phase shift is satisfactorily reproduced by the 
discretized solutions for continuum states. 
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In §2, we explain the formalism for treating the level density in the complex 
scaling method. We study the reliability of this method by applying it to a simple 
potential model in §3, and we investigate 4 He+n and a+a systems in §4. In §5, a 
summary and conclusions are given. 

§2. Continuum level density in the complex scaling method 

Here we briefly explain the complex scaling method (CSM). In the CSM, the 
spatial coordinate r and the wave number k transform as 

17(0) :r^r exp(t0), k -> k exp(-z0), (2-1) 

where U (0) is a scaling operator and is a real number called a scaling parameter. 
Under this transformation, the asymptotic divergent behavior, ~ exp(i/c r r), of a res- 
onant state with a complex wave number k r = k — i*y is changed into a damping form, 
exp {i(n — ry)(r cos0 + ir sin 0)} = exp {(7 cos0 — Ksin0)r}-exp {^(kcos + 7sin0)r} 
for > tan -1 j/k. Therefore, resonant states and bound states are obtained as dis- 
crete solutions of the complex scaled Schrodinger equation 

H(9)$ e = E$ e , (2-2) 

where H{6) = U(6)HU~ 1 (6). Because we require the complex scaled Hamiltonian 
to have no singularity, the scaling parameter has an upper limit, 9c- For the 
Gaussian potential, 9c = vr/4. For 6 < 9c-, the solutions of bound states and 
resonances with 7/re < tan0 are square-integrable, because of their damping forms 
in the asymptotic region. Therefore, employing an appropriate scaling parameter 9, 
we can derive resonant states in addition to bound states using a square-integrable 
basis expansion, for example, in terms of harmonic oscillator or Gaussian functions 

N 

<2> e = ^c n (0)<^. (2-3) 

71=1 

In Fig. 1, we present a schematic eigenvalue distribution for the complex scaled 
Schrodinger equation. It is seen that the energies of bound states are not changed 
from the spectral positions of the original Hamiltonian. The eigenvalues of the res- 
onant states, which are particularly noteworthy, are obtained as E = E r — ir r /2, 
where E r and r r are the energy and width of a resonance, respectively. By contrast, 
the continuum spectra of the Hamiltonian H{6) are distributed on the 20-lines orig- 
inating from every threshold. If we do not apply the complex scaling, the original 
Schrodinger equation gives the continuum spectra, including resonances on the pos- 
itive energy axis. Under complex scaling, the resonances for which 7//C < tan0 are 
separated from the continuum, and the rotated continuum spectra starting from dif- 
ferent threshold energies are separately obtained on different 20-lines. Furthermore, 
when we apply a basis function method to solve the complex scaled Schrodinger 
equation, these continuum spectra are discretized on different 20-lines, as shown in 
Fig. 1. 



4 



R. Suzuki, T. Myo and K. Kato 



Im(Energy) thresholds 





\ Re(Energy) 


K » • *-J- 






V 29 / V29 / 








%^ ^-^^^ 












X. v— -i' resonances 




•t resonances *l 




rotated continuum states 



Fig. 1. Schematic energy eigenvalue distribution for a complex scaled Hamiltonian. 



Let us return to the problem of the level density. The level density p(E) of the 
Hamiltonian H is defined as 



p{E) =£6{E - Ei), 



(2-4) 



where the quantities E\ are the eigenvalues of H, and the summation and integration 
are taken for discrete and continuous eigenvalues, respectively. This definition of the 
level density can also be expressed using the Green function: 



p(E) = --Im<! Tr 

7T ' 



1 



E-H 



= Im I dr ( r 

IT 



E-H 



(2-5) 



Here, applying the CSM and the extended completeness relation (ECR) 8 )' 9 ) to the 
expression of the Green function, we obtain 



p(E) 



-— Im J dr (^r 
-— Im J dr (^tq 

-— Im / < 

vr J 



i 



i 



E-H{9) 



E-H 
re 



u{e)- 1 u{6) 



y ^(r)<%(r) 
^ E-E B ^ 



R 



E — E R J Lg 



E — Ei. 



(2-6) 



where Nb and are the numbers of bound states and resonances in the wedge 
region between the real energy axis and the 20-lines, respectively. A detailed expla- 
nation of the extended completeness relation is given in Ref. 8). 
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In the integration over r in Eq. (|2-6|) . the bound state and resonance parts 
are easily found to be unity, because of the normalization of the wave functions, 
but the continuum part cannot be calculated, due to the singular nature of the 
integration. This singularity is eliminated when we discretize the continuum spectra 
using the basis function method with a finite number N of basis functions. Then, 
the approximate density of states p$ (E) for the basis number N is expressed as 

N B n r 

p ^E)=J2S(E-E B )--ImJ2 w -- 

B R 

N-Nb-N^ 

-> E w^wy < 2 ' 7 > 

As explained above, the energy of the resonance is obtained as Er = E r — ir r /2, 
and thus each resonance term has the Breit-Wigner form 

i -r r /2 , s 

Im = — . (2-8 

E — Er (E-E r f + r?/4: y > 

For the continuum part, discretized continuum states are obtained on the 20-line 
in the complex energy plane, = £j? — i£[, where £{/£ k R = tan 29. Therefore, 

the continuum part in the level density can be expressed in terms of a Lorentzian 
function whose form is similar to the Breit-Wigner form: 

1 -£ I 

Im — = * T . (2-9) 

E-£ k {6) {E _ £ RY +£ ^ 1 1 

Inserting Eqs. (|2-8j) and (|2-9j) into Eq. (|2-7j) . we obtain the level density in the basis 
function method as 



Vr; , - V ;, j., 2 j N-N B -N« R £j 

n^(E-E r ]i + r r y4 + K E {E _ £ R)2 + s i 



2 • 

(2-10) 



Here, it is noted that (E) has a dependence, but p{E) does not. This 
dependence problem of p^ (E) is due to the fact that we employ a finite number of 
basis functions, and it can be solved by introducing A{E) defined in Eq. (|lTj) . The 
continuum level density (CLD) A(E) is expressed as a balance between the density 
of states p(E) obtained from the Hamiltonian H and the density of continuum states, 
po(E), obtained from the asymptotic Hamiltonian Hq in the form 

A(E) = p(E)- P0 (E), (2-11) 

where p(E) is defined through subtraction of the bound state term from p(E). Phys- 
ically, A(E) represents the density of unbound levels, which result from the inter- 
action with a finite range. This can also be understood from the fact that A{E) is 
related to the phase shift caused by the interaction. 
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In the basis function method with a finite number N of basis states, we have 

A»{E)=p${E)-p» {e) {E). (2-12) 

The first term on the right-hand side represents the level density in which the bound 
state term is subtracted from Eq. (|2-7jl . and the second term is expressed in terms 
of the eigenvalues £2(9) = £2 R ~ of the asymptotic Hamiltonian Hq(9), which 
has only continuum spectra on the 2#-lines: 

i N cOI 
p N r E ) = ±y tk (2-13) 

Thus, we have 



K r /9 N-N B -Ni 

nA$(E) = Y- Lr i l 9/ + y 



( S _ + r 2 /4 Z_. {E -£*)* + £* 



v c07 



-V ^ (2-14) 

V(^-^r) 2 +^ oj 

As shown by the numerical results presented in the next section, the 9 dependence 
of Ag (E) disappears through the cancellation of the 9 dependence in the second 
and third terms of Eq. (|2-14jl . When we consider a small value of 9, and therefore 
no resonance exists in the wedge region, the CLD can be expressed in terms of only 
the second and third terms. 

§3. Simple potential model 

We now examine the reliability of the present method for a simple potential 
model. As a schematic potential, we employ the CGKPM potential, 13 ) whose res- 
onance structure has been studied in detail. The Hamiltonian in this case is given 
by 

H = T + V : T = -— V 2 , F(r) = -8.0e- ai6r2 +4.0e- a04r2 , (3-1) 

where we set h 2 / fi=l (MeV/fm 2 ) for simplicity. The Schrodinger equation for this 
Hamiltonian is solved by applying the basis function method, and thus we write 

N 

i>(r) = ]T Rt(r)Y em (r), R,{r) = £ c%(r, b t ). (3-2) 

Ira i 

For each partial wave, we use Gaussian functions 14 ) with different size parameters 
as basis functions: 



(r,bi) = N e (bi) ■ /exp 



2bf . 



2 e+2 



^^-'^WTTWT^ • (3 ' 3) 



1/2 
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Fig. 2. Energy eigenvalue distribution of the 
1~ states for the complex scaled Hamilto- 
nian of the simple potential model given by 
Eq. IjHTfl . The circles represent eigenvalues 
and dashed lines are 2#-lines. 
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where the parameters {foj : i = 1, 2, ■ ■ ■ , N} are given by a geometrical progression 14 ) 
of the form 

bi = b Y- 1 . (34) 

Here, &o and 7 are the first term and the common ratio, respectively. We employ 
N = 30, 60 = 0.2 fm and 7 = 1.2 in the following calculations. Of course, the same 
results are obtained even if other kinds of basis functions (for example, harmonic 
oscillator functions) are used. 

In Fig. we plot the energy eigenvalue distribution of the 1~ states. One bound 
state (—0.67 MeV) and many resonant states exist: The three lowest resonances are 
1.1710 - i0.0049 (MeV), 2.0175 - i0.4863 (MeV) and 2.5588 - il.7378 (MeV). 15 ) 
The lowest resonance is obtained with the CSM for 8 = 5°, but the second lowest 
resonance is not obtained for this scaling parameter value. The second resonance 
appears clearly when 9 > 10°. The continuum solutions vary slightly from the 20- 
line, and the dispersion increases for large values of 9. However, this distribution of 
continuum eigenvalues depends on the choice of the basis functions, and it creates 
no difficulty in the CLD calculations. 

Using these eigenvalues, we calculate the level density pff (E) given by Eq. (|2-7|) . 
and we display the result in Fig- EI Oscillatory behavior is seen at 9 = 5°, but this 
oscillation is smoothed when 9 is larger than 10°. Even at 9 = 5°, the oscillation 
may disappear if we employ a large number of basis functions so that the intervals 
between the discretized continuum eigenvalues become smaller than their imaginary 
parts. However, it is easier to choose a larger value of 9 in order to increase the 
imaginary parts of the discretized continuum eigenvalues. The intervals between the 
discretized continuum eigenvalues depend on the number N of basis functions. The 
critical value of 9 may be defined as the scaling angle at which the imaginary parts 
of the discretized continuum eigenvalues become larger than the intervals between 
the eigenvalues. This critical value of 9 depends on N, and therefore we express it 
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Fig. 4. The complex scaled CLD Ag(E) for 9 = 10°, 15° and 20° and the exact solution obtained 
by solving the scattering problem. 

as On- When becomes larger than 10° in the present simple potential case, p^ (E) 
exhibits the same behavior, and therefore we can set 9n ~ 10°. For 9 > 9n, only 
the absolute values of pff (E) depend on 0. 

This dependence of the absolute values of pQ (E) can be canceled through 
subtraction of p^ e (E); that is, we show that the CLD (E) defined in Eq. (|2-14|) 
has no 9 dependence for > n- In Fig. |3J we plot the CLD (E) calculated 
for = 10°, 15° and 20° and compare it with the result of the exact calculation. 
Here, "exact" means that we calculate the CLD A(E) from the phase shift using 
Eq. . The phase shift is obtained with the help of the scattering solution without 
any approximation. From Fig. |1J we see that it is quite difficult to distinguish the 
plots of A^ (E) calculated for 9 = 10°, 15° and 20°. They are all consistent with the 
exact calculation. This result indicates that the CLD A(E) can be approximated 
by A^ (E) in the CSM, and the phase shift can be obtained from A^ (E) without 
solving the scattering problem. 



We now apply the present method to realistic two-body systems of 4 He+n and 
a + a. The 4 He+n system has rather broad resonances but no bound state, and the 
a + a system also has no bound state but a sharp resonance due to the Coulomb 
barrier. The Coulomb potential is a typical long-range potential and is represented 
by the asymptotic term of Hq. The antisymmetrization among clusters in both the 
systems 4 He+n and a + a is carried out with the orthogonality condition model 
(OCM). 16 ) We show that the present method is very useful in analyses of continuum 



§4. Applications to He+n and a + a systems 
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states of such realistic cluster systems. 
4.1. 4 He+n system 

The wave function of 5 He with spin J is expressed in the 4 He+n cluster model 

as 

$ J ( 5 Ke)=A{<l>( 4 Re)-tf el (r)}, (4-1) 

where .4, <?( 4 He) and V^i( r ) are the antisymmetrizer, the internal wave function of 
4 He assuming a (Osi/2) 4 configuration, and the relative wave function between He 
and the valence neutron, respectively. We solve the relative wave function ^^{t) by 
applying the OCM. This yields 

[T rel + V an {r) + A |0pf>(0pf| - E] ^ J cl (r) = 0, (4-2) 

where T re i and V an (r) are the kinetic energy and potential operators for the 4 He- 
n relative motion, respectively. In this calculation, we use the so-called KKNN 
potential 17 ) for V an (r), which provides an accurate description of the low-energy 
scattering data for this system. The third term, constituting the non-local potential 
in Eq. (|4-2|) represents the projection operator to remove the Pauli forbidden (PF) 
states [which is the (Osi/2) state in this case] from the relative motion, 18 ) and A is 
taken as 10 6 MeV in this calculation. 

Equation (|4-2j) is solved by using the basis functions, as explained in the previous 
section, and we obtain 

N 

^(r) = [Ye(f)xi/2]jW(r), <pt{r) =J24<t>l(r,bi), (4-3) 

i 

where [YHr)xi/2] j ^ s a f unc ti° n of the orbital angular momentum and spin coupled 
to J, and the radial wave function <pe(r) is expanded in the Gaussian basis functions 
{(j>e{r,bi)}, which are defined in Eq. (|3-3|) . 
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Fig. 5. Energy eigenvalue distributions of the 4 He-n system for the J n = 3/2~,l/2~ and 1/2" 1 
states, where 6 is taken as 35°. 
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Fig. 6. Continuum level densities of the 4 He-n system for the J" = 3/2 , 1/2 and l/2 + states. 



Using the same basis set as in the case of the simple model, we calculate the 
energy eigenvalues of the complex scaled Hamiltonian with 9 = 35°, and the results 
for the three states 3/2~, l/2~ and l/2 + are shown in Fig. [5J We can see that each 
of the states 3/2~ and l/2~ has one resonance pole, corresponding to the observed 
resonances of 5 He. The 1/2+ state has no resonance. Resonant structures of 5 He 
have been investigated in detail with the complex scaling method by Aoyama et al. 19 ) 
In addition to resonances, the discretized continuum solutions have been obtained 
along the 2#-line. Several continuum solutions are off the 20-line. It is believed that 
the reason for this is that the couplings between the continuum and resonance are 
not correctly described because the number of basis functions is not large enough. 
However, the resonant solutions are obtained with appropriate accuracy, and the 
CLD is obtained from these continuum solutions satisfactorily, although the positions 
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Fig. 7. Scattering phase shifts of the 4 He-n system for the J n — 3/2~, l/2~ and l/2 + states. 
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of some continuum solutions are slightly off the 2#-line. 

Applying Eq. (|2-14j) to the obtained eigenvalue distribution of the complex scaled 
Hamiltonian for the 3/2 _ , 1/2 - and l/2 + states, we calculate the CLD of the 4 He- 
n system. The results are shown in Fig. |HJ It is seen that each of the 3/2~ and 
1/2 - states has a peak, but the l/2 + state has no peak. The position and width of 
the peaks in the CLD for the 3/2~ and l/2~ states corresponds to their resonance 
energy and width. These results are very similar to those for the CLD distributions 
calculated by Arai and Kruppa, 4 ) except for the absolute strengths. In the calculation 
of the CLD carried out by Arai and Kruppa, the result depends on the smoothing 
parameter. 

To see the reliability of the CLD obtained here, we calculate the phase shift from 
the obtained CLD. In Fig. [7J we show the scattering phase shifts of the 3/2~, 1/2 - 
and l/2 + states. We compare these results with the exact phase shifts, and we find 
very good quantitative agreement between them for every state. 

4.2. a + a system 

Similarly to the above, we now calculate the CLD and the scattering phase shifts 
of the a + a system. The important point in the calculation of the a + a system 
is the treatment of the Coulomb interaction. Because the Coulomb interaction has 
a long-range nature, we must include the Coulomb interaction in the asymptotic 
Hamiltonian Hq. 

The relative motion between two a clusters is described within the OCM as 

V r J el (r) = 0, (4-4) 

where V^ a and V^ a are the folding Coulomb and nuclear potentials obtained by 
assuming a (Os^) 4 harmonic oscillator wave function with oscillator constant v a {= 
4f£) = 0.2675 fm" 2 for an a cluster, respectively. When we employ the Schmid- 



Trei + V£(r) + V£ a (r) + A £ |<Apf)(0pf| - E 

PF 
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Fig. 9. Scattering phase shifts of the a-a system for the J n — + , 2 + and 4 + states. 
Wildermuth force 20 ' as the two-nucleon force, they are expressed as 




Vbexp 



u a /j, 2 
r 



v a + 3^/4 



(4-5) 
(4-6) 



where erf(x) is the error function, and Xjj = 2.445, Vo = —72.98 MeV and // = 
0.46 fm~ 2 are the folding parameter, the strength, and the range parameter of the 
Schmid- Wildermuth force, respectively. The fourth term in Eq. IJ4-4JI is the projection 
operator to remove the Pauli forbidden states (the 05", IS" and 0D states in this case) 
from the relative motion, 18 ) and A is taken as 10 6 MeV as well. We solve the complex 
scaled Schrodinger equation Eq. (|4-4[) in the same way as we solved the Schrodinger 
equation for the simple potential and 4 He-n systems. Using the obtained eigenvalues 
for J w = + , 2 + and 4 + , we calculate the CLD. In the a-a system, however, the 
eigenvalues of the asymptotic Hamiltonian Hq must be solved with the Coulomb 
potential: 

4e 2 

H = T rel + — . (4-7) 
r 

The results of the CLD are shown in Fig. |SJ They have a sharp peak corresponding 
to the resonance in each state. This result is quite similar to the results by Arai and 
Kruppa, 4 ) in which case the smoothing was performed in terms of several smoothing 
parameters. Better agreement is obtained for a narrower smoothing parameter. 

Integrating the obtained CLD, we obtain the scattering phase shifts. The results 
are shown in Fig. |SJ The scattering phase shifts are nearly identical to those obtained 
from the scattering solutions. The resonance width of the + state is very small in 
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comparison to the resonance energy. For such a case, it is necessary to carefully 
integrate the CLD to obtain accurate phase shifts. These results indicate that the 
present method to calculate the CLD is also very powerful even for a long-range 
interaction, such as the Coulomb potential. 

§5. Summary and conclusion 

We have shown that the level density is properly described in the CSM with a 
basis function method. In the expression we obtained for the level density, the ex- 
tended completeness relation of the CSM plays an important role, and it divides the 
level density into three terms, i.e., bound states, resonances and continuum states. 
We investigated the approximate description of the continuum states in terms of dis- 
cretized eigenstates that are obtained through diagonalization with a finite number 
of basis functions. Furthermore, it is not necessary to use a smoothing technique, 
such as the Strutinsky procedure employed by Kruppa and Arai for the singular 
level density arising from the discretization of continuum states. In the CSM, con- 
tinuum states are expressed in terms of eigenstates of complex eigenvalues along the 
rotated branch cut with the angle 26, and the Green function of the continuum part 
is expressed as a sum of Lorentzian functions rather than delta functions. There- 
fore, no singularity appears. This result indicates that the CSM provides a very 
powerful method for discretizing continuum states. The discretization of a contin- 
uous function using Lorentzian functions would be understood through comparison 
to the wavelets 21 ) that have recently been developed as a powerful tool facilitating 
transformations between analogue and digital data in information science. 

The level density smoothed in the CSM, however, has a dependence on the 
scaling angle 6, because a finite number of basis functions is used in the approximate 
description of the continuum states. We showed that the continuum level density 
(CLD) in the CSM, obtaining by subtracting the level density for the asymptotic 
Hamiltonian, is independent of the scaling angle and consistent with the exact CLD. 
These results indicate that we can calculate scattering phase shifts or S-matrices from 
the CLD obtained by solving an eigenvalue problem in the CSM with a finite number 
of basis functions. We found that this method is quite effective in the treatment of 
4 He+n and a + a systems without and with the Coulomb interaction, respectively, 
which were previously studied by Arai and Kruppa. 4 -* 

Considering the successful results of this method for simple two-body systems, 
it would be interesting to apply it to coupled-channel systems and three-body sys- 
tems. For coupled-channel problems, the extended completeness relation providing 
the foundation of this method has been proven in the framework of the CSM. 10 ) 
Therefore, it is conjectured that the present method will be effective here too. How- 
ever, the three-body problem is still open. 
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